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A MATRIC GENERAL SOLUTION OF LINEAR DIFFERENCE 
EQUATIONS WITH CONSTANT COEFFICIENTS 


Ralph W. Pfouts and C, E. Ferguson 


Linear difference equations with constant coefficients have been used 
for many years in various disciplines, and a general method of solution has 
become well known. In this paper an alternative general solution is pre- 
sented; a statement of the linear difference equations system in terms of 
linear vector spaces then follows naturally from the alternative solution. 
The special case of the second-order linear difference equation is consid- 
ered first, not merely because it is a simple case but also because it is a 
form which has many applications in the social sciences, especially in 
economics and psychology. 

1. The Second-Order Linear Difference Equation.—The homogenous 
portion of the second-order linear difference equation 


(1) Y,=@,Y F 
may be shown as 

Yt 
(2) 


0 OILY 
Equation (2) clearly shows the recursive nature of a difference equation 
system. However, all of the essential features of the system are retained 


if (2) is abbreviated to read 


Yt 


(3) = 


Using R to designate the square matrix in (3), we note that the char- 
acteristic equation of F is 
(4) |R-Al| = A?-a,A-a,=0, 
which is also the reduced characteristic equation of the square matrix in 
(2). It is obvious that the right-hand member of (4) is also the auxiliary 
equation of 

Since the trace of a matrix equals the sum of the characteristic roots, 


and since the determinant of a matrix is the product of the characteristic 
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roots, we may rewrite (3) as 


Yt “A 

(6) 
0 

Taking advantage of the recursive nature of (3), we write the general homo- 


geneous solution as 


1% 
0 
where y, and y, are the initial values of the dependent variable expressed 
as deviations from the value of the particular integral. Then if ) is a par- 
ticular integral of (1), a general solution of the equation is 


(8) == 
1 0 

where Y , and Y, are the initial conditions. Y, is shown as a function of ¢ 
in (8). 

A question of importance in applications is whether (1) is stable, in 
the sense of economics, or damped, in the sense of physics, That is, we 
wish to ask whether Y, converges to - It is clear from (8) that the solu- 
tion is stable if, and only if, R converges to the null matrix. This, in turn, 
can occur if, and only if, the A’s are less than unity in absolute value, a 
point which is proved in the next section. 

Frequently it is desirable to determine stability from the coefficients 
of (1). In the second-order case, the necessary and sufficient conditions 
are 

1-|a,|>0, 


(9) l-a,-a,>0, 
l+a,-a,>0. 


These conditions are well known, cf. for example [6]. 
2. The General Linear Difference Equation. — The matric general solu- 
tion of the general linear difference equation 


(10) 


may be found in a manner analogous to that used for the second-order equa- 
tion. Thus the homogeneous portion of (10) may be exhibited as 


| i 
lye 
| 
| 
| 
| 
| 
| 
| 
| 
n 
| Y t a; t—-2 
] 
| 
| 
| 


t—1 


L Ye—n+ 0 0] LYn—t J 


Using R to designate the square matrix in (11), it follows immediately 
that 


If one sets the expression of the right in (12) equal to zero, it is clearly 
the auxiliary equation of (10). Let us call this expression equation (12a). 
From the elementary theory of equations, we recall that the coefficient a; 
is equal to the sum of the combinations of the n roots of (12a) taken i ata 
time. Thus we write 


(13) S; 


to indicate that ¢ roots are included in each product and that there are (“) 


terms in the summation. 
Then the solution of the homogeneous system (11) is 


0 

(14) : =} : 


where / is the unit matrix and y,, --+, y,_, are the initial values of the 
homogeneous variable. Similarly, if ' # is a particular integral of (10), the 
general solution of (10) is 
] t—nt+1 


Let us turn again to the question of stability. We first prove a necessary 
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and sufficient condition which depends upon the characteristic roots of FR. 
Next we present one necessary and three sufficient conditions; these con- 
ditions are developed directly from the elements of FR. Finally we state the 
so-called Schur-Cohn conditions which are both necessary and sufficient 
for stability. These conditions also involve only the elements of k. How- 
ever numerical values of the elements must be known if the Schur-Cohn 
conditions are to be useful. 

Theorem: A necessary and sufficient condition for stability of the 
difference system (14) is that the characteristic roots of R be less than 
unity in absolute value. 

The proof of this theorem is well known cf., for example Mirsky, [4, 
pp. 328-29]. But we can sketch a simpler proof using the information at 
hand.’ In equations (3) and (14) the full difference systems are abbreviated 
so that the minimum and characteristic polynomials of R are equal. Hence 
R is nonderogatory and, for this reason, similar to the Jordan canonical 
orm 
(16) J=A+C, 
where A is an (nxn) diagonal matrix whose elements are the characteristic 
roots of R and C is an (nxn) matrix with ones on the first superdiagonal 
and zeros elsewhere. Consequently R* converges to the null matrix if J°* 
does so. Finally, if R* converges to the null matrix as s increases without 
bound, the solution of the difference system converges to the particular 
integral. 

Consider the binomial expansion of J* for s > 7: 

For any integer which equals or exceeds n, the power series expansion 
contains exactly n terms since C” is the null matrix. Without loss of gen- 
erality we may assume that A = A/. Substituting this expression in (17) 
yields 


To prove sufficiency, we assume |A| < 1. As s increases without 
bound, A°+0; thus the main diagonal of J* contains only zeros. The re- 
maining terms in (18) give the entries on the successive superdiagonals of 
J*, namely + ae for the &-th superdiagonal. Each of these terms goes 
to zero as s increases without bound. This is seen by considering the gen- 
eral term (ener as a sequence of numbers generated by increasing the 
value of s. It is well understood that this sequence converges to zero as s 
increases without bound. Hence all other elements of J* approach zero, 
and J* converges to the null matrix. 

Necessity is immediately proved by assuming |A| > 1. Then |A|*>1 
and the first term in (18) does not vanish. Accordingly J* does not converge 
to the null matrix. This completes the proof for the Jordan canonical form. 
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Since F is similar to J, it also proves the same property for R*. 

Necessary conditions for stability involving the coefficients of (12) 
may be obtained from (13) by noting that if the A’s are of modulus less than 
unity, each term in the summation in (13) will be less than unity, and con- 
sequently 


(19) (7). 1,2 


The inequalities of (19) are necessary conditions for stability. 

Sufficient conditions for stability have been developed by several dif- 
ferent writers. Some of these conditions may be easily shown. The left- 
hand side of the determinental equation in (12a) may be written as a func- 
tion of A: 


n n—1 

The zeros, or roots, of this polynomial occur, of course, at the points at 
which the polynomial is itself equal to zero. 

As early as 1799 Gauss [2] showed that the polynomial in (20) has no 
zeros outside of the circle |A| = r, where 

r= max (1, /2 S$), 

and S is the sum of the positive C; coefficients. For stability the roots 
must lie within the unit circle. Consequently, the Gauss test is not di- 
rectly relevant since the minimum possible value of r is unity. That is, 
one or more roots could lie on the unit circle, giving rise neither to ex- 
plosive movements nor damping. However, the Gauss test does show that 
the sum of the positive coefficients in the characteristic equation must be 


less than 2~* for stability. 
A somewhat similar, but more stringent, test was constructed by Walsh 


[8]. He showed that the roots of (20) lie within the circle given by 


n 


In particular, stability is assured if 


LE 
j=l 


This provides a rapid test based upon all of the coefficients of (20). 
A more easily applied test was developed by Kojima [3]. He proved 


that all of the roots of (20) lie within the circle |A| < r, where 


r= max 


for j = 1, 2,+++,2-1. Using this test, the sufficient conditions for stability 
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may be quickly determined by scanning the coefficients to see if either of 
the coefficient ratios exceeds unity. A similar test was suggested by Mon- 
tel [5], but it does not discriminate as finely. 

We have now shown some necessary and some sufficient conditions, 
There exists one set of necessary and sufficient conditions which rely 
only upon the values of the elements of &. These conditions were originally 
developed by Schur [7]; they were later refined by Cohn and others, and 
since that time, the conditions have usually been called the Schur-Cohn 
conditions. Specifically, both writers proved that it is necessary and suf- 
ficient for stability that the matrix 


1 


ol 


. 
— 


= 1, 2, +++, n) 
be associated with a constrained positive definite quadratic form. 

3. The Difference Equation System in Linear Vector Space. The matric 
general solution of the linear difference equation system suggests, in a 
natural way, a statement of the linear difference equation in terms of linear 
vector spaces. To approach this subject we define 
Yi 7 


Viens 1° 


where the V’s are column vectors belonging to the polynomial in (12), Thus 
each V, is an n-dimensional vector in the homogeneous variable of (10). 
Making use of the matrix R, we may show a sequence of these vectors as 


(23) V, RV, R7V,, 
The vectors of (23) are merely a consecutive sequence of n of the vectors 
defined in (22). 

Any consecutive sequence of n vectors of (22) or (23) are linearly in- 
dependent. This may be demonstrated by noting that the elements of any 


| 
a 0 eee a 
n n—1 1 
1 
a eee a eee a 
} 1 n N—t+2 
a: a: eee 0 «62 
(21 D,= 
a eee a a; 
n 1 t—1 
a a ] ms 
ie 
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vector in the sequence are related by the homogeneous difference equation 


and the elements of any consecutive set of n vectors are similarly related. 
Consequently scalar multiplication of any one vector in a consecutive se- 
quence means that all other vectors in the sequence are multiplied by the 
same scalar. Thus it is impossible to induce linear dependence in (23). 
From the definition of a homogeneous difference system, we observe that 


Thus the (n+1)-st vector of a consecutive sequence is a linear combination 
of the preceding vectors. Consequently a set of (n+1) consecutive vectors 
is linearly dependent, but any set of n consecutive vectors is linearly in- 
dependent, 

The arguments just presented demonstrate that n consecutive vectors 
of (22) or (23) may be considered as a basis of a linear vector space I" of 
dimension n, This space we call the solution space. 

Now apply the theorem which states that any vector in the set span- 
ning I’ may be expressed uniquely as a linear combination of vectors con- 
stituting a basis of I’. Hence any solution vector V, in I’ is expressible 
uniquely as a linear combination of consecutive vectors. 

We now digress to prove the following 
Lemma: If W,, W,, +++, W, denote the n linearly independent column vec- 
tors of a non-singular n-square matrix P, if A is an n-square matrix, and if 
AW = b tons then the column vector 5, by 
bn jl is the j-th column vector of P~'AP. 

By hypothesis the n vectors of P are linearly independent. Thus the 
vector AW ,, which is the j-th column vector of the matrix AP, may be ex- 
pressed as stated in the lemma. If we denote the i-th row of P~' by W*, we 


have 


WW = 


where 6,_; ; is the Kronecker delta with respect to ¢ and j. Then it follows 


immediately that 


; = 6 


t~i,j* 
This proves the lemma [1]. 

We return now to the chief argument. By definition of a homogeneous 
linear difference system 
(25) V,=RV,_, =n, «+, t) 


From (24) and (25) we obtain 
(26) RV, 


Now by virtue of the lemma proved above and (25) and (26), we see that a 


; 
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vector V; may be taken as a column of an n-square matrix P such that 


P~'RP = R’= 


where ?’ is the transpose of FP. 
It is apparent that 


Yten+ 1 Yt-an+a2 


Thus a set of n consecutive vectors in |’ may be used in forming a con- 
tragredient transform of A, and this transformation yields the transpose 
of R. 

Rand RP’ are rational canonical matrices. We now introduce the theorem 
which states that if A is an n-square matrix and |A—A/| has only one in- 
variant factor different from unity, the matrix is similar to a rational ca- 
nonical form, And under these conditions the invariant factor in question 
is the characteristic function of the rational canonical matrix[1, pp. 170-1]. 

As a consequence of the theorem stated in the preceding paragraph, 
the right member of (12a) 1s the only invariant factor of |R-—A/| which is 
different from unity. Since (12a) is the auxiliary equation of the linear dif- 
ference system, we may also say that the auxiliary equation is the only 
nontrivial invariant factor of the system. The auxiliary equation of a linear 
difference system is usually treated as if it were merely a convenient trick 
for obtaining the general solution. The arguments of this paper show that 
it is much more than this. Indeed it is fundamental to the difference equation 
system, for it is the only invariant factor different from unity and it also 
contains the constants associated with a basis of the linear vector solution 
space, 


FOOTNOTES 


The authors are indebted to Professor Leonard Carlitz of Duke Univer- 
sity for suggesting this method of proof. 


The argument of this section relies on that of [1] pp. 207-11. 
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ON UNEQUAL PARTITIONS OF INTEGERS 


Iwao Sugai 


SECTION I. INTRODUCTION 


The ancient Greeks attached considerable religious and philosophical 
significance to numbers. Among all, Pythagoras highly valued square num- 
bers such as 4, 9, 16 ---, he expressed the square of any positive integer, 
say 7 by the sum of 7 gnomons which increase their magnitudes in the order 
of appearance.’ Thus, he partitioned i? into i distinct positive integers. 
This note is an extension of such an idea, i.e. for given positive integers 
say i and n, how many ways are there to express i” by i distinct number 
elements, possibly including non-gnomons. Thus the method taken in this 
note does not resort to the well-developed theorems on prime numbers and 
congruences, but it uses extensively arithmetical progressions and finite 
differences. The scope of the general parition of integers under such 
methods is so vast that two conditions are imposed to facilitate study. 
They are discussed in the following sections. 
NOTATION. 

For positive integers i and n; (i)” is expressed as the sum of i dis- 
tinct positive integers, such as 


m=t 
= Anim Where i = 1, 2,3, 1, «++ 
m=1 


From these number elements an infinite “right-angled-number-triangle” is 
formed. If the diagonal elements and horizontal elements, or if the diagonal 
elements and vertical elements have zero n-th finite difference and n > 3, 
the total possible number of such “number triangles” is at least 3-2"~°. 
Thus (i)” has at least 3-2"~* partition types under the aforementioned 
conditions. 


SECTION Il. 
CONDITIONS AND TWO INTERESTING TYPES OF PARTITIONS 


Condition 1 
For positive integers i and n, (i)” is expressed by the sum of i dis- 


tinct positive integers, i.e. 


m=t 
(1) ("= S im Where = 1, 2 3, 
m=1 


Eq. (1) is arranged as follows, 
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n,t,8 


where j = 7+ 1 
Condition 2 
In the diagram above the number elements must satisfy one of the fol- 
lowing conditions : 
(a) Diagonal elements An and horizontal elements 
-+-, 2) for all 7’s must have zero n-th finite differences. 
or 
(5) Diagonal elements A, ; ; and vertical elements A, i,p (p = 1, 2, «+, #) 
for all z’s must have zero n-th finite differences, 
In order to handle many varieties of the infinite “right-angled-number-tri- 
angle” (as it appears above), three basic classifications are defined at 
this point. 
TYPE (1) 
(I) It is generated by the expression 


m=t 
(i)” = [m” —(m—1)"] . 
m=] 


(II) Vertical elements A, ; p are the same for all 7’s if p is held constant. 
(III) Both diagonal and horizontal series are identical and they have n! for 
the (n—1) th finite differences. 
TYPE (2) 
(I) It is generated by the expression 
m=t 
(3) (= 14 
m=) 
(II) Number elements are in ever-increasing order, i.e. Ay, im 
any m, where j = 7+1. 
(III) Both diagonal and vertical elements have (n—1)! for the (n—1) th finite 
differences, 


‘ 
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TYPE (8) 
Any other varieties of partition types which do not belong to either 
Type (1) or Type (2). 
It is interesting to observe that the case of Pythagoras (n = 2) belongs 
to Type (1) and there are no other types. Thus this note deals with cases 
of n > 3 to study many types of partitions of i”. 


SECTION III. GNOMONS FOR CUBE NUMBERS 


The literature in the history of theory of numbers mentions often the 
gnomons for the square numbers, but there seems to be little mention of 
gnomons for cube numbers. Now gnomons for the cube numbers are illus- 
trated by the following diagrams. 


25-2345 


19415417419 


Figure 1 


The number triangle consisting of these “successive” gnomons is repre- 
sented below. 


| 
: 
39274941 
3+5 749411 
| | | | 
13415417419 
| 
i 
3 
2° =3+49 
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749411 
4° 13415417419 


5% = 214+234+25+4+27+29 


m=s+i—] 
> (2m—1) where s = 1+ 


2 


This was first reported in the literature by G. Polya.” Another type of rep- 
resentation by gnomons for the cube numbers is also shown here. 


| 
| 


| 
| 
| 
| 


Figure 2 
From these gnomons the number triangle 
« 147 
3° = 1474+19 
4° = 1+7+19+37 


53 1+74+194+37+61 


is formed, 

These represent two methods for dividing cube numbers into gnomons. 
As shown later, the “successive” gnomons are very useful in describing 
the Type (2) partition of this note. This idea of gnomons extends to any 
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higher powered numbers, so that i” is expressed by the sum of i gnomons, 
SECTION IV. 
CONSTRUCTION OF VARIOUS TYPES OF PARTITION 


(a) n=3 case 

The construction of three types of partitions is illustrated. Take i = 2 
first. 2* must be expressed by the sum of two unequal positive integers. 
There are three ways to meet this requirement, namely 1+7, 2+6, and 3+5. 
These three ways come from 


2° = 27.2 = (1+3)-(1+1) . 


lj + 3 3 


“ ” 
by the “successive 
gnomons 
3+5 2+ 6 
Next take 7 = 3, and notice that the ways the partitioning lines are drawn 
must be followed exactly. 


33 37.3 = (1434+5)-(14+141) 


+ 37+ 5 5 
1 + 5 lj+ 5 
by the “successive” , 3 45 i+ 314 5 
gnomons 
7+94+11 
A pattern similar to this sort of partition is shown in the book by Uspensky 
and Heaslet.® For all other i’s the patterns are identical. 
(b) n= 4 case 
The construction of six types of partitions are explained in like man- 
fl +7) 
or 
24 = 2°.2 = 42+ 6 (1+1) 


+19 3+ 9 +15 


ner, 


3 + 


3 +13 5 +11 
At this stage it is useful to re-examine the generating function for 


; 
or 
34+5 
| + 7 7 21+ 6 2 + | 6 
1+ 7 7 6 2 + 6 
a) 1 +15 2 +14 4 +12 6 +10 a 
5 3 +| 5 
| 
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Type (2). It has the lowest power index with respect to the running index 
m, and Type (2) gives 2* = 6+10, Therefore, any expression which gives 
2 = 7+9 does not involve m, i.e. by this generating function (i)" is not 
expressed by 7 distinct positive integers in general even though 2* is ex- 
pressed by two positive integers. From the view point of gnomons for the 
fourth powered numbers, the partition type of 24 = 7+9 is also rejected. 


Gnomons for the 2‘ are shown below: 


‘ 
24 = 


5 


| 


Figure 3 
22749 
— 
daa 
|. 
Figure 4 


But the last case violates the nature of gnomons. Gnomons increase their 
magnitudes in the order of appearance. 

Pythagoras’ original gnomons are arrays of ones with width of a unit 
and his gnomons do not skip any array of ones, but some of our gnomons 
shown for n = 3 and n = 4 cases have width larger than a unit and our gno- 
mons skip many arrays of ones. These may be called “generalized” gnomons. 


SECTION V. THE FINAL RESULT 


The result of this study is to show that the total number of unequal 
partitions of 7”, obeying Conditions (1) and (2) is at least 3.2"~%. 
For the case of n = k, k&: any positive integer, (1)* is partitioned into 


4 
; 
q | 
| 
| 
| 
| 
| 
: 
d 
x 
; 
i 
: 
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i distinct positive integers. The possible number for various partitions for 
i* (i > 3) is equal to the possible number for various types of partitions of 
2* into two unequal positive integers, since the patterns for iP (i > 3) are 


identical for the patterns for ?. 

Now 2" is partitioned into two positive integers from the highest power 
index of m, running index (Type 1) to the lowest power index of m, running 
index (Type 2). é 


(1)+(2*-1) — Type 1 


(2) +(2*-2) 


Collection of Type 3’s 


2* (3) +(2*-3) 


— Type 2 


From the last line, S is solved to be 3.2*-8, Thus, there are at least 
types of partitions. 

The auxiliary conditions for the n-th finite differences for diagonal 
plus horizontal or diagonal plus vertical number elements are obvious by 
proper application of calculus of finite differences‘ to various generating 
functions. 

Of all these 3.2"~* types of partition, the Type (1) has been applied 
for a practical purpose in extracting the n-th root of a positive integer on 


digital computers. ° 
SECTION VI. NUMERICAL EXAMPLES 
Due to limited space, only the numerical examples of cases n = 3 and 


n = 4 are shown. 
- 3 case 


m=t 


Type (1) —(m-1)*] 
m=) 


=] 


= 1+7 


= 27 =14+7+19 


2 G4 1474194 37 


6° = 216 = 1+7+194+37+614+91 


73> = 343 = 14+74+194 +614+914+127 


: 
. . 
: 
i 
< 
a 
1 
2°. 
2 
tk 
ay 
bd e 
‘ 
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m=t 


Type (3) i alm? -(m-1)"1 


1° =) 
2+6 
= 3+9+15 
4+ 12+ 20+ 28 
5+15+25+35+445 
= 6+ 18+30+ 42+ 54+ 66 


7421+354+49+63+77+91 


27-345 

33 749411 

13415417419 

55 214+234+25+27+29 

6° = 314+334+354+374+39+41 

73 = 434454 47+ 49451453455 


. 
. 


m=t 
(1) i*= [m4 —(m-1)"] 


=l=1 
=14+15 
1+15+65 
= 14+15+65+175 
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7) 
= 
m=t 
Type (2) = 
1 3 _ l 
n= 4 case 
Type 
Lund 
4 
94 
34 
44 
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54 = 625 = 14+154+65+ 175+ 369 


m=t 
Type (3)-a = ~(m—1)*] 


m=1 
24=2+14 
34 = 3421457 
4* = 44+28+76+ 148 
54 = 54+354+95+ 285+ 305 


m=t 
Type (3)-—b = + (3m?-3m+1)] 


m=1 
1f=1 
24 = 34+13 
34 = 74+25+49 
= 13+43+79+121 
54 = 214+67+119+177+241 


m=t 
Type (3)—c i‘ = > [(2m—1)i?] 
m=1 


1*=1 

24 = 4412 

34 = 94+27+45 

4* = 16+48+80+112 

54 = 254+754+1254+175+225 


3 
. . 
; 
2 
4 
. . 
. . 
. . 
| . 
. 


is 
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m=t 
Type (3)—d_ = (i(é+1)(2m—1) ~ (3m ?~3m+1)] 
m=1 


=] 
5411 

= 11+29+4] 
19+53+81+ 103 
29 +83 + 131+173 + 209 


m=t 
Type (2) i* = ili? -i+1+42(m-1)] 
=] 
1*=1 
24 = 6+10 
34 = 214+27+: 
44 = 52+60+68+ 76 


54 = 1054+1154+1254+1354+145 
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ROCKERS AND ROLLERS 


Gerson B. Robison 


Some years ago, while picking up my small son’s toy blocks, I became 
intrigued with the possibility of finding a cylindrical surface upon which 
a plank would roll in neutral equilibrium. The results appeared as a prob- 
lem in the American Mathematical Monthly. (v. 59, Oct. 1952, p.551) The 
solution is repeated below. 

For some reason, the normal mathematical impulse to generalize lay 
dormant until recently, when it was stimulated by a proposed science ex- 
hibit at New Paltz Teachers’ College. There is a straightforward approach 
to the problem which yields several interesting and constructible models. 
In this approach, the upper, moving surface, (the roller) is defined in polar 
co-ordinates, with the center of gravity at the pole, while the fixed sur- 
face, (the rocker), is defined by rectangular co-ordinates with the y-axis 
positive downward. 

The requirement of neutral equilibrium means, first, that the center of 
gravity of the roller must travel in a horizontal path and, second, it must 
remain directly above the point of contact of the two curves in all posi- 
tions. In addition, the roller must actually roll into each position. These 
conditions suggest the following formulation : 

In Fig. 1 the roller curve is given by r = rs), 0 = 0(s). The pole is at 
0’, and the parameter s is the length of arc A’P. The rocker has for its 
z-axis the line of travel of the pole of the roller. Its curve is given by 2 = 
a(s), y = y(s), s here being the corresponding arc length AP, equal to A’P. 
We assume that all four functions are differentiable. In the figure, the ini- 
tial ray and the origin have been chosen so that 2(0) = 6(0) = 0, but this is 
actually a boundary condition which may be modified in some cases, 


Rocker RPS 
Roller 

Initial ray of roller 
S° Common tangent TPU 


Figure 1 
139 
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The condition of neutral equilibrium now gives us 
(1) ns) = y(s) for all values of s 
and the geometry yields 


dz r dé 
(2) —=t PU = . 
d 


The special points where = a = 0 do not give trouble in general. 
8 8 
To illustrate the method, let us start with the roller whose surface 


cross section is the straight line r = sec 6. By differentiation we get od = 
r 


1/n/r?-1. Using (1) and (2) above we obtain =1/Vy?-1. Sor =cosh™ 'y+C. 


We now take 0(0) = 2(0) = 0, which amounts to a correlating of the two co- 
ordinate systems. Then r(0) = y(0) = 1, and hence C = 0. We thus arrive at 
the equation y = coshz for the rocker. (Fig. 2) 


Figure 2 


One can start with a rocker equation, although this seems to be less 
frequently fruitful of interesting results. Consider the family of oblique 


lines: y = kx. For these rockers = 1, whence = 1, leading to = 
dy dr ft 


Inr+C. In this case it is convenient to take @(0) = 0, but : for 2(0). Then 


y(0) = (0) = 1, again making C = 0. The roller family is then given by r = 
eX 9, (Fig 3) We have tabulated below some of the more interesting pairs. 
In each case the c.g. of the roller is at the pole of its equation. The last 
three columns correlate the two curves. 

Number 5 has some special interest, for the two parabolas are actually 


. 
i 
2 
5 
f 
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congruent. In fact, it can be shown by synthetic methods, using the focus- 
directrix definition of the parabola, that it is the only curve for which the 
rocker and roller can be congruent and in contact at corresponding points. 


Figure 3 


Referring to fig. 1, if F 1s the point of the rocker corresponding to 0’ of 
the roller, then FP = PO’, and P is therefore on the parabola determined 
by 00’ and F, 

In actual construction of the models, some practical problems arise. 
Obviously, it is desirable to have surfaces with a high coefficient of fric- 


tion, since the value of — on the roller must not exceed this constant, in 


operation. A layer of fine sandpaper or of felt is effective, but allowance 
must be made for its thickness in the layout. 

A thornier problem is that of getting the center of gravity where it be- 
longs. The plank of number 1 is easy, and the parabola of number 5 will 
have its center of gravity properly located if it is a segment of altitude 
1 2/3. However, number 9 calls for a circular roller with the center of grav- 
ity on its circumference. One way of dealing with the problem is to double 
the roller. A system of two tangent circles has the centroid properly placed. 

In number 8 of the table, we find that a cardioid rolls on acycloid. 
However, both curves have cusps, which means that there will be not only 
excessive slope, but also interference in the rolling action. By jumping 
from one lobe of the cardiod to the other along the common tangent, we 
obtain a roller which will roll on a combination of cycloid and catenary, 
the join taking place at the points of 60° slope. In a similar fashion, a 
catenary can be combined with the circular rocker of number 9 at the 45° 


Catenary arc — PS 


Figure 4 
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TABLE OF ROCKERS AND ROLLERS 


ROLLER 


Equation : r = Description 


Straight line 


equiangular spirals 


lt+e 


(e> 1) hyperbolas 
1+ecos 


l+ecos 


2 


parabola 
1+cos 


spiral 


je hyperbolic spiral 


1-kcos 6 limacon (cardioid if k= 1) 


2 cos 6 circle, unit radius 


r= cosné “roses” 
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ROCKER 


Equation: y 


Description 


coshez 


catenary 


oblique lines 


catenary-like 


—ecos \/ ——< 


l-e 


l+e 


sinusoidal 


paraboia 


(2>0 for k>0) 
((14+ ka) (2>0 for -1 <k<0) 
(2 <0 for k<-1) 


parabola-like 
hyperbola-like 
parabola-like* 


exponential 


z=t—-ksint 


y= 1-kcost 


trochoid (cycloid) 


=4 


circle, radius 2* 


yy? = 


ellipses* 


ber 6. 


* Rocker concave upward in number 9, number 10, and third species of num- 
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points, for a roller consisting of two tangent circles and their common ex- 
ternal tangents. (Fig. 4) 

We conclude with a comment on parametric forms. If 6, 2, and y ( = r) 
are differentiable functions of any parameter ¢, then (1) can be reduced to 
dz 
dt 
may be used to obtain 2(¢), y(¢) being given with r(¢). In particular, we may 
take @ as the parameter, whereupon the foregoing equation becomes simply 


“* = r, Number 8 is an illustration of this. Conversely, given parametric 


dt 
equations for the rocker, to obtain parametric representation of the roller 


we may use 29 _ 1 


dt y dat 


= <. Hence, if the roller is given parametrically, the above equation 
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TEACHING OF MATHEMATICS 
Edited by 


Joseph Seidlin and C. N. Shuster 


This department is devoted to the teaching of mathematics. Thus articles on method- 
ology, exposition, curriculum, tests and measurements, and any other topic related to 
teaching, are invited. Papers on any subject in which you, @8 a teacher, are interested, 
or questions which you would like others to discuss, should be sent to Joseph Seidlin, 
Alfred University, Alfred, New York. 


USE OF MATRICES IN TEACHING 
CONIC SECTIONS 
A. R. Amir-Moéz 


Algebra as a tool has been used in problems of geometry. But we are 
always far behind in using a new tool when discovered. There is little 
hesitation in using determinants in elementary algebra. I don’t see why we 
should hesitate to use matrices in freshmen analytic geometry. I believe it 
brings a new interest and also it eliminates some tedious arithmetic. 

1. Definition: We are used to writing a point of the plane as P = (2, y). 
Let us ignore “ ,” and write it as (2 y). We call (2 y) a row matrix corres- 


x 


ponding to the point P. To the point P also corresponds (*'), called a col- 
umn matrix. Now we generalize this idea as follows: 


An array of numbers, for example 


3 2 
1 0 


0 


is called a matrix of two rows and three columns. In what follows we do 
not use matrices having more than three rows or columns. 
It is custommary to write a matrix as 


Two subscripts are used to indicate to which row and which column an 
element belongs. For example a,, is the element located in the second 
row and in the third column. 

Two matrices are equal if and only if all corresponding elements are 
equal, 

2. Addition of matrices: We define addition only for matrices which 
have the same number of rows and columns, The example 
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23 


: 

: 

: 
a a a 
ne 11 12 13 

a a a ah 

21 22 23/ - 

3 


shows how addition is defined, i.e., corresponding elements have been 
added. 

3. Multiplication of a number and a matrix: To multiply a number p by 
a matrix we multiply every element by p, as in the following example 


a a 
Pay, 


Pa. 


4. Product of two matrices : Two matrices can be multiplied if the num- 
ber of columns in the first one is the same as the number of rows in the 
second. We shall demonstrate the multiplication by the following example 


=| a. b +a..b a a_.b 
211 

b 

1 


22 1 22 21 1919+ 22 22 


° 1 82 21 31 12 32 223 


3 


1 
1 


b 


21 


is obtained by adding the product of the first element of the third row of 
the first matrix ard the first element of the first column of the second ma- 
trix to the product of the second element of the third row of the first matrix 
and the second element of the first column of the second matrix. 

It is important to observe that the operation of taking the product of 
two matrices is not necessarily commutative. For example 


2 1\/-l -2 -5 -6 -l -2 21 -6 


0 o} \-4 o/,\-3 \o 2 -3 }. 
For convenience we shall denote a matrix by a capital letter such as 


$). 

For matrices A, B, and C the reader may verify that A(B+C) = AB+AC 
if the multiplication and addition are possible. 

5. Transpose of a matrix: If we interchange rows and columns of a 
matrix A we call the new matrix the transpose of A and we denote it by A’. 


For example if 


A=/4,, @,,\, then A’ = %,} 
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a a 
res 19 11 12 
a a c c 
81 32 31 82 
We see that, for example, oe 
81 
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It is easy to show that (AB)’ = B’A’, (note the order). We leave it to the 
reader to verify this property. We also observe that (A+B)’= A’+ B’. 

6. Straight line: An equation of a straight line which passes through 
two points (2 ,, y,) and (z,, y,) is 


y-y, = ,). 


This equation can be written as 


Y~-Y 


where ¢ is the common value of these fractions. Sometimes it is desirable 
to write these equations as 


= 2,-2,) = tl 


(1) 


ty ,) = tm 


which are called parametric equations of the line and 7 and m are called 
direction numbers of the line. Now we write the matrix equation of the line 
which is essentially the same as (1); 


x 


(2 y)=(2, y,)+t m). 


Now let X = (@ y), X, = (x, y,), and D =(1 m). Then (1) can be written as 


X=X,+tD, or X’=X%+0D". 


7. Conic sections: An equation of the second degree in two variables 
is a conic section and can be written in the form 


ax? +2bry+cy? +2pr+2qytr = 0. 


This equation can be written using matrix multiplication in the form 


(1) (x y b \f 

be }= (0). 
The reader may multiply and check the equation. Now let X =(z y 1), 


Q=/a bp 
b c/, and beg 


3 
1 
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Then (1) can be written as 
XQX’=0. 


8. Intersection of a line and a conic: Let XQX’ = 0 be a conic section 
and let X = X ,+¢D be a straight line. Here X = X ,+¢D is taken to be 


(2 y y, 1)+tl m 0). 


The points of intersection are the solutions of the system of equations 
XQX’°=0 
X=X,+. 

Setting X = X ,+¢D in the first equation we get 

(1) (X ,+tD)Q(X 5 +tD’) = 0, or X ,QX QD’+ tDQX + t?DQD’=0. 


(We have to keep in mind that multiplication is not commutative). We leave 
it to the reader to show that 


X QD’ = 
Thus (1) may be written as 
(2) (DQD’)t? +2(X QXi=0. 


We see that (DQD’), (X,QD’), and X ,QX{ are numbers (matrices of one 
element). We also note that to a point we can make correspond either (z y) 
or (z y 1). We shall denote these row matrices by X and D = (1 m) or 
(2 m 0). We shall use them where suitable. 

Let us study the following : 


DQD’=(1 m O)fa b = (7 
be 


m 


bec 
p q tf 


Thus (2) can be written as 


(3) (DAD’)t? + ,QD’)t+X 


Now the values of ¢ obtained from (3) can be substituted into X = X ,+¢D 
and we get the points of intersection. 

9. Discussion of the intersection : 

I. If DAD’ # 0, then the line intersects the conic in two points, real 
or complex conjugate, or coincident. This is determined through the dis- 
criminant of [8.(3)]. 

II. If DAD’ =0, X ,QD’ #0, the line and conic intersect in one point. 


Ill. If DAD’ = 0, X ,QD’=0, but X ,QX 40, the line does not intersect 


: 
5 
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the conic, 
IV. If[8.(3)] is an identity, then the line is actually a part of the conic. 


10. Center of a conic: A point (z,, Yo) is called a center of a conic if 
for any point M of the conic, there is another point N on the conic such 
that P is the midpoint of the segment MN. 

Now to find the center let P = (2, y, 1) and M =(@ y 1). The equa- 
tion of the line through P and M is 


(1) X=P+tD, 


where D = (x-2, y-y, 0) [see 6.]. As in [8.] the points of intersection 
of (1) and the conic XQX’ = 0 are found by 


(2) (DAD’)t? + APQD’)t+ PQP’=0. 


Let DAD’ # 0; otherwise there will not be two points of intersection. If 
t, and ¢, are the two roots of (2), then the two points of intersection will be 


M = (2, y) = (a) +¢ yy + ¢ 


and 


Thus 


N =(2,,y,) = (2) +t, yy + 


t itt, 


2 


). 


P = (25, Yo) = +[a-2,] 


t +t 
This implies that __? = 0. But t +t, is the sum of the roots of (2) and we 
have 
QPD" 
DAD’ 
Since DAD’-4 0 the equation PQD’ = 0 gives the centers if there are any. 
That is 


(a 


be qghy-y, }=9, 


Pq 0 


or 
(ax,+ by, + = 0. 


Now two cases may be considered : 
I. If (x, y) = (@,, y,), then every point of the conic is a center and the 


conic is a double line; 
Il. If (w, y) # (2, y,), then we must have 


ba,t+cy,+qg=0, 
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¥, l),a 
cj=0. 
Pq 


This system of equations will give the centers, or we might write it as 


az+by+p =0 
(3) 
ba+cy+q=0. 
ll. Discussion of the center: The pair of equations [10.(3)] can be 
considered as a system giving the point of intersection of two lines. 


I. If b?~ac = 0 and ag-pb # 0, then the two lines of [10.(3)] are par- 
allel and they do not intersect. Thus the conic does not have a center 
(parabola has no center). 


II. If 6? -ac= 0 and ag— bp = 0, then the two equations of [10.(3)] are 
the same and there are infinitely many centers (two parallel lines). 


If b?-ac #0, then there is a unique center, 
12. Examples: 
(i) Find the center of 


52x? -72ry + 73y?+82-294y-1167 = 0. 


Since (52)(73) -(~36)? # 0, a center exists. The equation can be writ- 
ten as 
(x y -36 
-~36 73 -147Ny)=0. 
4 -147 -1167/ \1 
The center is the solution of 
(x y 1); 52 -36 
70 0, 
4 -147 
522 —36y+4 = 0 
—362+ 73y-147 = 0. 


The center is (2, 3). 
(ii) Find the center of 


1627+ 242y + 9y?+ 1502 -200y- 1000 = 0. 


Since (12)?—(16)(9) = 0, there is either no center or infinitely many cen- 
ters. We have 


(x y 1)/16 12 75\ fe 
2 9 -100)y 
75 ~-100 -1000 


150 
or 
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But p/q # 16/12, i.e., ag-pb # 0. Therefore there is no center. 
(iii) Find the center of 


64a? + 36y? + 4802 - 360y+675 = 0. 


We see that (64)(36) —-(-48)? = 0 and a/b = -64/48 = p/g = - 240/180, i.e., 
aqg-pb=0. Therefore there are infinitely many centers. We expect the conic 
to degenerate into two lines parallel to 


642 -48y +240 = 0, or = 0. 
We shall get them as follows. The two lines must be 
8x-6y+h=0 and 8r-6y+k=0. 
Writing 
(8x-6y+h)(8a-6y+k) = 642? -96ry + 36y? + 4802 -360y+675 =0, 

we obtain 

h+k = 60 

hk = 675, 


and fh and & are the roots of x? ~602+675 = 0 which are / = 45 and & = 15. 
Thus the conic is 


(82-6y+45)(82-6y+15) = 0. 


13. Change of coordinate system to the most convenient position: If 
a conic section has a center we translate the origin to the center. This 
translation changes the conic into the form 


az? +2bay+cy? =1 
because the origin becomes the center of symmetry. Now by a suitable ro- 
tation we change this to the form Ar?+ py? = 1. This is sometimes called 


the canonical form. 
14. Change to canonical form: Let a conic be of the form 


(1) ax*+2bay+cy?=1, 


Then (1) can be written as 


Now we know that the equations of a rotation of the coordinate system are 


= 2 cos 0-y,sind 


# sind 


This in matrix form is 


-sin@ cos@/. 
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Let 88) and x =(@ y) and X, = (2 


-sin@ cosé@ y ,)» Then we see that 


1 


sin 6\;/cos @ 0 


)-1- 8. 
01 


-sin@ cos @/\sin@ cosé 
Therefore 
sin@ cosé@ 01 


We also see that /A = A] = A, where A is any matrix. Now (1) can be writ- 
ten as 


(2) (7) = XAX’= XIAIX’ = XS’SAS’SX’= X SAS’X*.. 
Now the problem becomes that of finding S such that 
°) =B. 
Op 


This problem happens to have a solution, and we are going to solve it. We 
see that SAS’ = B can be written as SAS’S = BS, or SA = BS, or 


-sing@ cos Alt i cos 6/. 


This equation gives 
acos 0+ bsin@ = Acos@ -—asin@+bcos@ 
bcos -bsin@+ccos@=cosé. 
These equations are written as 
(a—d) cos 6+ bsind = 0 
cos 6+(c-A) sind =0, 
6-bcosé =0 


Since sin@ and cos @ can not be zero at the same time we must have 


(a—A)(c-A)-b2 =0, or =0, 


nd 
(a~p)(c-p) - 62 = 0, or 


We observe that A and » must satisfy the same equation, i.e., 
(4) b6| =0, or =0. 
c-2 


This equation always has real roots because its discriminant is 
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(a+c)?-4(b?—ac) = (a—c)?+4b7>0. 


We also see that A # y; otherwise a = c and } = 0 which implies that (1) is 
a circle. Obtaining A and » we substitute in (2) and we get 


X SAS’X (2, =(1), or Ae? + py? =l, 
0 

Indeed there are two rotations which change (1) to canonical form. Thus 

the other possibility is pr?+Ay? = 1. 


Now we did not need to find the rotations of coordinate system. But in 
case the rotations are desired we use equations (3) and find the rotations. 


15. Example: Let 2277+ 42y+5y? = 1 be the conic for which the coor- 
dinate system is to be rotated in order to eliminate the ry term. Clearly 


is the same conic section. Thus [(14.(4)] will be 


and we get A = 1, p = 6. We have either 
v?+6y?=1, or +y?=1. 
Now suppose we want to find the rotations. Equations [14.(3)] will be . 
(2-1) cos 0+2 sin@ = 0 
cos 6+(5-1)sin@ = 0, 
(2-6) 6 = 0 
sin 0+(5-6)cos@=0. 


These equations are all equivalent since A and ,» satisfy [14.(4)]. We see 
that any of these equations reduces to 
cos @+2sind=0, or tan@=-1/2. 
Since we are interested in sin@ and cos @ we observe thatcos? @= ——__ = 
l+tan76 


4/5, and since < 0 < 2/2 we have cos = 2/5/5 and sin@ = - 5/5. 
Therefore the rotation which changes the conic into 2?+6y? = 1 is 
2/5/5/, 
= 
y = (-/5/5)e, +(2V5/5)y, 


or 


3 
| 
: 
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Now for the other rotation we write A = 6 and y» = 1. It is sufficient to use 
only 


(2-6) cos 6+2sin0=0, or -2cos@+sing=0, or =2. 


Comparing this with the previous case we see that one rotation takes the z- 
axis into a direction perpendicular to the direction into which the other 


rotation takes the z-axis. 


1 ‘ 
Now cos? @ = = 1/5 and as in the previous case, - 7/2 < 
l+tan? 6 


7/2 and we have 


cos6=\/5/5 and sin@ = 2/5/5. 


Thus the rotation which changes the conic into 627+ y? = 1 is 


-2/5/5 /5/5/. 
16. Example: (Reduction of a conic with center); By a translation and 
rotation simplify 


522? = 0. 
The reader may see that the center exists. We write the equation as 
(x y  ~36 Ay = 0. 
~36 73 
4-147 -1167’ \i 
The center is obtained from 
4 (2 y 1l)/52 -36, =0, 
-36 73 
4-147 
52¢-36y+4 = 0 
+ 73y-147 = 0. 


: i The solution is (2, 3). Therefore 


9 
(+2 


y=y,+ 3 
is the translation which puts the origin at the center. We observe that the 
substitution for « and y does not change the form of the second degree 


terms and the first degree terms will disappear. We have to compute only 
the new constant term. This constant is found as follows: 


52(2 7) — 72(2)(3) + 73(3 2) + 8(2) — 294(3) - 1167 = - 1600. 


= 
ey 
| 
i 
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Since the left side of the equation equals the constant term when 2, =y, =0, 
5207-722 y,+73y? = 1600. 
Now we proceed as in [15.]. 


52-2 -36 
-36 73-2 


gives \ = 100 and yp = 25, Thus either 
1002? +25y2 = 1600, or 2522+ 100y? = 1600, 


= 0 


which reduces to 
= 64, and +4y2 = 64 respectively. 


We leave it to the reader to find the rotations if desired. 
17. Reduction when no center exists: Let the conic be 


and b?-ac=0. 
Note that there might be infinitely many centers [see 11.II]. But we treat 
the problem in the same way. We write the equation as 
(x = (-2pr-2qy-r). 
be 


Now we use a rotation to change the left side to the form 


We know that A and ,» are the roots of x? —(a+c)r-(b?-ac) = 0 [see 14.(4)]. 

But using the relations between the roots and coefficients of this equation, 
we see that A+p = a+c and Ap = 56?-ac = 0. Thus either A = 0 or p = 0. 
Suppose » = 0 and dA # 0. This is possible since » = 0 and A = 0 imply that 
a+ec = 0 and thus b?-ac = b? +a? = 0 which means a = 0, b = 0 and it is of 
no interest. Therefore we have 


(1) 


A=at+ec and p=0 


and (1) will be 
(a, 5) z, =([fe,,y,)]. 
0 0 


In order to determine f(z,,y,) we have to find the equations of rotation. 
We use [14.(3)], i.e., 
1 a® 


asin@—bcos@=0, or tan@=6/a. Then cos? @ = = 


39 43,32 
But 5?~ac = 0. Thus l+tan*@ a+b 


2 
a a 
cos” 6 = 
2 a+c 
a“+ac 


cos 0 = 


a+c 


Since -7/2 < 6 < m/2 we have 
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We observe that 5? = ac > 0 implies that a and c have the same sign and 


this shows that a/(a+c) > 0. Now siné = Subs 6 or sind = as £ 


Let us study different cases: . plana’ 
I. Let 6 > 0. Then 4 = and sin = 


(i) If a > 0, then we have sin @ = pve 
(ii) Ifa < 0, then 


(ac)(a) 
ate a*(a+e) 


sin @ = 


a 


(ii) If a > 0, then we have sin@ = -y oe 


(ii) If <0, then we have sind The result is 


—(-a) a*(a+c) a+c 
cos 


and the sign for the latter equation is the same as the sign of ab. 
18. Example: Simplify 


+ + 9y? + 1502 1000 
by a rotation and translation of coordinates. 


We see that cos = sinfd=y 3 Thus the rotation 
9+16 94+16 5 


4/52,-3/5y, 
3/92 ,+4/5y,, 
and A= 16+9 = 25, » = 0. Thus we have 
252? + 150(4/5 ,-3/5 y ,) ~200(3/5 2 ,+4/5y ,)- 1000 = 0, 


or 25a,?-250y,-1000=0, or 2? = 10(y,+4). 
By a translation of coordinates we get 


2 
= 10y,. 


Queens College, New York 
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POLYNOMIALS AND FUNCTIONS 


Louis O. Kattsoff 


1. Introduction : 
The purpose of this paper is to clarify the’ nature of a variable and a 

function. This will be done by starting with the definition of a polynomial 

over a set, and defining the concept of a variable and function in terms of 

procedural rules, 

2. Polynomials : 

Let us consider a set S containing elements a; (¢ = 1, 2, 3, ---) and op- 


erations ‘o’ and‘{]’, We take any symbol ‘z’ and any othersymbol ‘+’ which 
are to be used as follows: 

2.1. The symbol z will be juxtaposed to any element a; on the right of 
a,, to form an expression we will call a monomial. 

2.2 If x is repeated after a given a, then the number of x’s in the mon- 
omial may be indicated by writing a superscript to z. 
Example 2.3. @,2, @,, etc. are monomials. 


wx may be written as a,,2°. 


5 


5 


3 


a, may be writtenas a,,2°. 


The sub-scripts on the elements of S are solely for the purpose of distin- 
guishing the elements of S, 
2.4. The symbol ‘+’ occurs only between monomials. For example: 


8 57 
907% 


2.5. Definition of a polynomial over a set S. A polynomial in a symbol 
2 over a set S, is an expression of the form 


0 1 2 3 n 


where Ay € S. 
3. Comments : 3.1. 2° is usually not indicated. 

3.2. Since z is not necessarily an element of S the right hand jux- 
taposition of z to an element of S cannot be taken to indicate the multipli- 
cation of the element of S by the z. a,x" is not an element in S, even if S 
happens to be closed. 

3.3. The + is neither the o nor the [] of S. Therefore the polynomial 
in 2 over S cannot be identified with some element of S even if S happens 
to be closed, 

3.4, We can form a polynomial in a, over S, a, « S, as follows 
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2 8 n 
2g + +@,a, t+ eee + 


a,+a,a 
and this will not be an element in S either. 
4. Operations with polynomials in 2 over S. 
Definition 4.1. A polynomial in z over a set S of elements {s,, Sorts 


8, +++} will be denoted by S(z). 


3 n 
S(z) = 8) +8 + 


Definition 4.2. If S (x) and S,(z) are polynomials in z over a set S 
such that 


2 
S (a) = Sot +82 


= 


then S (x) = S,(z) if and only if sx" = six", 

Definition 4.3. If S(z) is a polynomial in z over a set S, then the de- 
gree of S(x) will be the largest n for which x” has a coefficient s, 4 0. 

Definition 4,4, The power of a monomial is the exponent of z. 

It follows that S (x) = S,(z) implies s; = s/, and the degree of S (2) = 
degree of S, (2). 

Definition 4.5. If =a, and +5 2+ 


b,x? +-..+b,2" are polynomials in x over a set S with operation [J (addi- 
tion) anda;«S, S, then 


S (x) @S,(2) = (a, (a,b, +--+ (a,b, 


Theorem 4.6. If S (x) and S,(x) are polynomials in 2 over a set S with 
operation [7] defined in S, then S,(z) @S,(z) is also a polynomial in z over S, 
Proof: By definition 


S (x) @S,(2z) = (a,(5,)+(@ + (a, 
But since a; and 6; « S (and we assume S closed with respect to (7) 
Let 


Therefore @S,(z) = c, 
5. The generalized Polynomial: Since nothing is said about the nature of 
the indeterminate z, we can form polynomials in y over a set S with oper- 
ation |], also. However, we would not necessarily want a polynomial in « 
over S to be distinct from every polynomial in y over S. To justify treating 
two polynomials as the same we cannot use a definition but will need to 
make use either of a procedural rule, or a rule of replacement. 

Rule of Replacement 1. If S(z) is a polynomial in z over a set S with 
operation |7] defined, then if the symbol 2 is replaced everywhere it occurs 


by another symbol of the same kind, for example, y, the resultant polynomial 


+C 
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S(y) will be equivalent to S(z). E.g. 
54274298 +24 = 54 
The same result may be achieved by 
Procedural Rule 1. If S(z) is a polynomial in z over the set S and 


S(y) is a polynomial in y over the set S, with the corresponding coeffi- 
cients of S(z) and S(y) equal, then we will proceed as if S(z) = S(y). The 


symbols x and y must be of the same kind, 
E.g. ¢ because ‘2’ and‘7’are not sym- 
bols of the same kind. If the two polynomials are taken over the set F of 


real numbers FR, butz¢ R. 
The Rule of replacement evidently follows from the Procedural Rule. 
Procedural Rules will be directives to behave in a certain fashion. 
The essence of the concept of function lies in these procedural rules, as 


we will show. 
Theorem 5,1. If S,(z) and S (y) are polynomials over S in x and y with 


corresponding coefficients equal. 
S (a) = S,(y). 


Proof: S ,(y) differs from S (x) only in the fact that whenever there is 
an z in § (2), there is a y in S,(y). Therefore S,(z) and S,(y) have the same 
degrees and the same corresponding coefficients. By Procedural Rule 1, 
we are entitled to assert S (z) = S,(y). 

This suggests that we need to generalize the notion of the indeter- 
minate so that various types of symbols can be used to replace it. 

Definition 5.2. We adopt a meaningless symbol ¢. A polynomial in ¢ 
over a set S with the operation |7] defined in S is an expression of the form 


0 


Procedural Rule 2. In a polynomial S(2) over a set S, when the symbol 
¢ is replaced by any other symbol, the resultant expression will also be 
treated as a polynomial over S, 

Definition 5.3. If P(Q) and @() are two polynomials over a set S_ with 
operation defined in S, then P(2) = Q() if and only if for every symbol w, 
P(w) = Q(w). 

6. More comments. 

If Cis replaced by x, we obtain a polynomial in 2 over the set S. 

If Cis replaced by 2, we obtain a polynomial in V2 over the set S. 

We have used |_| to denote the operation of “addition” in S. We esta- 
blish 

Procedural Rule 3. In a polynomial in x2 over a set S with operation 
(|, replace the |]by + whenever it occurs. 

Since S (x) and S,(z) are not polynomials but names for polynomials, 


“S (2) +S$,(2)” is the name for the sum of the two polynomials. 
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Evidently the + between two monomials in a polynomial is not the same 
as the + which replaces the [] in S. No confusion results however from the 
indiscriminate use of these two symbols. 

7. Mappings. 

Evidently for any replacement of ¢ by some symbol w in S(), we obtaina 
mapping @: w + S(w), where S(2) is a general polynomial over the set S 
with operation +. 

Definition 7.1. If S(Q is a polynomial in ¢ over a set S with operation 
+ defined in S$, then S(Z) will be called a polynomial function over S. 

A polynomial function over a set S defines a set of polynomials over S. 

Procedural Rule 4. To obtain polynomials over a set S, we replace the 
meaningless symbol in S(Z) by indeterminates. 

Suppose we have a polynomial in z over a set S. Given S(x)=a, +a ,2+ 


+a,27+ +a,2". 
We consider what happens when we replace z by elements of S. Let 


s ,¢S. Then by our procedural rule we can write S(s,)=a,+@,8,+@,8°++--+a,8/. 


S(s,)is not yet an element of S because although a; « S and s, € 8, 
s, ¢ S since juxtaposition was not equated with multiplication in S. 

If S is a set with only one operation (e.g. a group), then S(s) cannot 
be an element of S. 

If S is a set with two operations, addition and multiplication, then we 
can adopt a 

Procedural Rule 5. If S(z) is a polynomial in z over a set S with addi- 
tion and multiplication in S defined, then we can identify the + between 
monomials with addition in S and juxtaposition with multiplication in S, 
whenever z in S(z) is replaced by an element s; of S. 

On the basis of P.R.5, we can establish the following mapping 


Since we did not assume S to be closed S(s,) is not necessarily an ele- 


ment in S after the operations are performed. 
If we add the condition that S be closed, then we have 


+ S(s;) where 


8; + S(s;) > 8; where s,, 8; € S. 


It follows that 6 maps S into itself. 

We can now define a function: 

Definition 7.2. A function is a polynomial in z over a closed set S 
(with addition and multiplication defined in S), such that only elements of 
S are used to replace z; or S is mapped into itself by the polynomial when 
z is replaced by elements of S. 

Definition 7,3. A variable in a polynomial in z over a set S (with addition 
and multiplication defined in S) is an indeterminate in the polynomial 
which may be replaced in accordance with procedural rules, by elements of S. 
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MISCELLANEOUS NOTES 
Edited by 


Charles K. Robbins 


Articles intended for this department should be sent to Charles K. Robbins, Depart- 
ment of Mathematics, Purdue University, Lafayette, Ind, 


RECEDING GALAXIES 
Lindley J. Burton 


Hubble proposed a generally accepted astronomical theory that gal- 
axies recede from each other at a rate proportional to the distance between 
them. The red shift of the spectral lines supports this theory. 

Here I wish to show with the help of vectors that if A and B each re- 
cede from 0 at a rate proportional to the distances 0A and OB, then B re- 
cedes from A at a rate proportional to AB in Euclidean space. 


Let 0A = and 0) = and 0B = and v(B,0) = tr. Then 
o(B, A) = resulting vector from ~tr, and tr, at B. Clearly AAOB ~ ABCD 
since <0= <C and the sides including these angles are proportional. Hence 
|BD| = t|AB\, and <CBD = <0AB, so that BD is an extension of AB. Thus 
BD = tAB, which proves the result with the same constant of proportion- 
ality. 

This result means that there is no unique center from which the gal- 
axies recede. Suppose that the centers of the galaxies are at the vertices 
of a grid infinite in extent and that the edge of the basic square of the grid 
is increasing in length at a fixed rate. Then the galaxies recede from each 
other at a rate proportional to their distance, 
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A NOTE ON SPACE CURVES 
Major H. S. Subba Rao 


Problem No. 248 appearing in Vol. 29, No. 1, September 1955, this 
Journal, proposed by Huseyin Demir, Turkey, can be extended to three di- 
mensions, 

It is well known that if |’, 
one correspondence such that the tangents at corresponding points are par- 
allel then the principal normals and binormals at corresponding points are 
also parallel. The two curves are said to be deducible from each other by 
a Combesure transformation. 

It can be stated as a theorem that the lines joining corresponding 
points on two twisted curves deducible from each other by a combesure 
transformation generate a developable surface and that the line joining the 
centres of curvature at corresponding points intersects the generator at a 
point on the edge of regression of the developable surface. 

Proof: Let fs = f(s ,) and *, = o(s,) be the vector equations to the two 


and Cy are two twisted curves in one-to-. 


curves, r, and r, being the vectors to corresponding points, s, and s, the 


arc lengths measured along the curves. If t, and t, are the unit tangents at 
corresponding points then =t,. If n,, are 


1 
the unit principal normals at the same points we have —! = —! and —2 - 


_* where p, and p, are the radii of curvature of the two curves. Since — = 
Po 
—— showing that Nn and —— = P 1/P o> Using 
Py ds, ds, 
Serret-Frenet formulae it will be obvious that the binormals are also par- 
allel. 
Let d be the unit vector in the direction of the vector r,—r,, that is, 


along the line joining the corresponding points on the two curves. If d’ = 


(a), b= we get by differentiating the vector equation 
s 
1 


ie 
4 
‘ 
wh 
id=r-r, (c= |r,—r,|) 
1 2 1 2 
— dr d 8 d 
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= 
= (ln), since ¢, =¢ 
Py 


1 . 


Squaring both sides and noting that d-d’ = 0 we get 


Differentiating the equation 


p 
OF (1-2) cos 6 = 
Py 


where @ is the angle between the vector d and a 
Hence we get 
Po \2sin? 


Py 


Similarly we can show that 
t,-d’=b= 


It is well known from differential geometry that the condition for the ruled 


surface generated by the vector r,-7, to be a developable surface is that 


b? a* sin? 


which is satisfied in this case. 
Writing the equation to the generator +, -r, as 


r =r, +kr,-7,) 


where K is a scalar variable, we get the point of intersection of this gen- 
erator with a neighbouring generator by solving for K from the equation 
dr dr, 


dr. r, ds 
+ K( 2 2 
de * de, 


=%,+K(t, -¢,£2) =0 
Py 


~ 
p 2 
2.2, 
1 
we get 
p 
(r -r,)-¢,1-—) = 
1 2 1 p 
1 
q 
(l1-—) 
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and substituting in the equation to the generator. The vector to the point 
of intersection C, which is a point on the edge of regression of the devel- 
opable surface, is given by 


r 


Since the principal normals at corresponding points are parallel, the 
centres of curvature and the generator lie on a plane. If M,, Mf, are corres- 


ponding points, C,, C, the centres of curvature at these points then it is 


easy to see that the lines MM, and C ,C, intersect externally at a point 


C’-such that 
C°M, : Py 


Thus the vector C’ is given by the same expression as for the vector to C. 


Hence C and (’’ coincide. 
It will be seen that the problem proposed by Huseyin Demir can be ob- 
tained by unrolling the developable surface into a p ane. 
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PROBLEMS AND QUESTIONS 
Edited by 


Robert E. Horton 


Readers of this department are invited to submit for solution problems believed to be 
new and subject matter questions that may arise in study, in research, or in extra-aca- 
demic situations. Proposals should be accompanied by solutions, when available, and by 
any information that will assist the editor. Ordinarily, problems in well-known textbooks 
should not be submitted, 

Solutions should be submitted on separate, signed sheets. Figures should be drawn 
in India Ink and twice the size desired for reproduction, 

Send all communications for this department to Robert E. Horton, Los Angeles City 
College, 855 North Vermont Ave., Los Angeles 29, California. 


PROPOSALS 


397. Proposed by Leo Moser, University of Alberta. 
Show that the Diophantine equation 


1 1 


has at least one solution for every n. 


398. Proposed by Huseyin Demir, Kandilli, Eregli, Kdz, Turkey. 


Determine the roots of the equations = 0, y? +2 =0 


where the coefficients (real numbers) in one equation are the roots of the 
other. 


399. Proposed by Nathan Altshiller Court, University of Oklahoma. 

Prove that the feet of the perpendiculars dropped upon the sides of a 
triangle from their respective Simson poles form a cevian triangle. That is, 
the lines joining those points to the respectively opposite vertices are 
concurrent. 


400. Proposed by Walter B. Carver, Cornell University. 
Given a point, a circle, and any curve in a plane. Construct an equi- 
lateral triangle having a vertex on each of them. 


401. Proposed by John M. Howell, Los Angeles City College. 

Given a sequence of numbers related by F(n) = aF(n-1)+5F(n-2), 
F(0) = c and F(1) = d, where n = 0,1, 2,-+- and a, 5, c, and d are any real 
numbers. Find a general form for F(n). 


402. Proposed by Chih-yi Wang, University of Minnesota. 
I. Show that if |2| <1, 
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F(z) = = 
n=k 


Il. Prove the same equality without using the binomial coefficients of 
negative arguments, 


403. Proposed by Vladimir F . lvanoff, Sar Francisco, California. 
Given four lines, a; (é=1, 2, 3, 4), determining a four space. Prove that 


if a; meets a;, , (mod 4), then all the given lines concur. 


SOLUTIONS 


Unorthodox Definitions 


376. [May 1959] Proposed by C.W. Trigg, Los Angeles City College. 

Identify the following unorthodox definitions: (1) the earthy or stony 
substance in which an ore or other mineral is bedded, (2) fully sufficient, 
(3) the four brightest components of @ Orionis, (4) type of speech fanci- 
fully exaggerated, (5) oddity, (6) a Franciscan friar, (7) adeposit formed 
in a liquid vegetable extract, (8) a large molding of convex profile, (9) 
fancied, (10) the series of air bubbles made by the breath from an otter 
under water, (11) a drain to carry off filthy water, (12) one sixteenth of a 
fluid drachm, (13) scheme, (14) depicted as walking, (15) absent-minded, 
(16) seekers, (17) complete, (18) invalid, (19) an interpreter of music. 

Solution by the proposer. (1) matrix, (2) abundant, (3) trapezium, (4) 
hyperbolic, (5) eccentricity, (6) Minor, (7) apothem, (8) torus, (9) ima- 
ginary, (10) chain, (11) sink, (12) minimum, (18) angle, (14) gradient, (15) 
abstract, (16) zetetics, (17) integral, (18) null, (19) exponent. Note the 
acrostic formed by the initial letters of the terms. 

Also solved by D. A. Breault, Sylvania Electric Products Inc., Wal- 
tham, Massachusetts. 


A Binomial Congruence 


377. [May 1959] Proposed by J. M. Gandhi, Thapar Polytechnic, Patiala, 
India. 
Prove that 


93n é 9 3n 
93 {mod 93 +3) 


i=0 
when n is of the form 2m, 4m+3 or 4m—1 respectively. 
Solution by Leonard Carlitz, Duke University. Put 


3n+] 
6N+2 yor 
( 


r=0 
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Then by the binomial theorem 


6N+ 2 


= (1+/3) +(1-/3 


2-3, so that 


=2, - 
and therefore 


8N+1, an+ 1 
2S = 2 (e ) 


Now if 7 is odd 


2 


r=0 
so that 
(1) § = (mod 2°"**) (n odd). 


If n is even 


2r<3n 


so that 


(2) S = 2(9n+1) (mod 2°"*4) (n even). 


Thus (1) and (2) furnish the desired result and indeed a little more. 
Also solved by Melvin Hochster and Jeff Cheeger (jointly) and the 
proposer. 


Sum of a Series 


378. [May 1959] Proposed by Barney Bissinger, Lebanon Valley College, 
Pennsylvania. 


Let a,, a, be arbitrary and n(n+1)a 


= n(n-1)a,-(n-2)a,_, for 


n+ 1 


n> 0. Find S = Xa, 


I. Solution by William M. Sanders, Mississippi Souther College. Since 


a, = a,/2 and ay = a,/6, a, depends only on a, for n > 2. By induction, 


a, = a,/n!. Consequently 


oo 


a,+a,+) —=a,+a, (e=1). 
n! 


16 
9 6N+2 
2S ) 
ete = 2+y3, 
9 , 
4 
| 
; 
é 
4 
S = a 
0 
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Il. Solution by J. L. Brown, Jr., Ordnance Research Laboratory, Penn- 
sylvania State University. From the given recurrence relation, we have the 
identity 


oo co oo 


n(n+l)a,, nin- la,Z" - (n= 12" 
n=] 


n= n=] 


which is equivalent to 


n=] 


Letting (Z) = 3 a2" , this latter identity may be written 
r= 


d? wt 3d KZ) 
2—, =z 40] -2 


d2 
dZ? aZ 


It is obvious that AZ + Be is the general solution of this equation, Fur- 
q 
d. 


ther, since (0) = a, an 


Z=0 


KZ) = (a,-a,)Z+a,e” 


The desired sum is equal to (1) = a, +(e-1)a,. 
Ill. Solution by Victor Ch’iu, Kent State University, Ohio. 


Mn+ Wan, , = 


n+ 1 1 


(n-1)na,, 


(n- 1)(n -2)a 


eee 


4-5a, = 4-3a,-2a, ‘ 


w 


3+ 2a,-a, 


2-3a. = 2-la,-0a, 


= 1-0a,-(-la, 


Me 
2 2 
— a = a - a 
dZ? nti dZ? n dZ 
ee: 
| 
ag 
a 
: 
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n(n+ lan, = -3a,-4a, 


n+1 1 


a,(2!)~" 


a,(3 


a,(4!)7! 


a,(n!)~! 


S=a 


a, 


0 2 3 4 


=a,+a,(e-1) 


If we assume a, = a,, then S = ea,. 


Also solved by George M. Bergman, Stuyvesant High School, New 
York; D. A. Breault, Sylvania Electric Products, Inc., Waltham, Massachu- 
setts; Leonard Carlitz, Duke University (two solutions); Jeff Cheeger, 
Erasmus Hall High School, Brooklyn, New York; Melvin Hochster, Stuyve- 
sant High School, New York; Philip M. Mills, Virginia Polytechnic Insti- 
tute; Paul Stygar, Norwich Free Academy, Connecticut; William B. Teach- 
out, Jr., Memphis, Tennessee; George B. Thomas, Jr., Concord, Massachu- 
setts; Chih-yi Wang, University of Minnesota; and the proposer. 


The Duel 


379. [May 1959] Proposed by J. M. Howell, Los Angeles City College. 

Two men fight a duel. They both fire at a given signal. If both are 
alive, they fire again at a given signal, and repeat the process until at 
least one of them is dead. If the probability that A kills B on any round 
is p, and the probability that B kills A on any round is r, find the proba- 
bility that A is alive, B is alive or neither is alive after: 1) n rounds, 
2) after an infinite number of rounds. 

Il. Solution by Glenn D. James, Los Angeles City College. If one 
knows that any one round is about to occur, the probabilities at the end 
of that round are: ; 


P(dd) = P(both dead) = pr 
P(aa) = P(both alive) = (1-p)(1-r) 
P(da) = P(b kills a) = r(1-p) 


1960) 169 
Brey 
Hence, 
2 
a,= 
3 
eee 
| 
Hence, 
=a,+a,(1+—+—+—+ 
21 3! 4! 5! n!} 
4 
3 
‘ 


MATHEMATICS MAGAZINE (Jan.-Feb, 


P(ad) = P(a kills 5) = pO =r) 


For convenience let & be the probability: P(dd)+P(da)+P(ad). Then the 
probability of death at the end of n trials is given on either side of the 
equation: 


kR = 1-P™aa) 


The desired probabilities are &P(ad), kP(da), and &P(dd), respectively. 
For n infinite, kis 1/R. 

Il. Solution by F. D. Parker, University of Alaska. Let P' = 1-p, r' = 
1-r, P”(AB) be the probability that both A and B are alive after the nth 
round, P’(AB') the probability that only A is alive after the nth round, etc, 
It is not hard to show by means of a tree or the appropriate stochastic 


1 1,1)” 1 1,1\7 
matrix that P"(AB) =(p'r')", P™( 4B) pu gipy 
l-p'r! l-p'r! 
1,1\7 1 1,1\7 
Then we have = P™B) = 
l-p'r l-p'r 
p'r+p(p'r')” prll—(p 'r*)”] 
l-p'r! l-p'r! 
p'r 
Il. Solution by the proposer. Consider: 


and P(A'!B') 


After an infinite number of 


rounds, these probabilities become , respectively, 


State 1 — Both alive 
State 2 — A only alive 
State 3 — B only alive 


State 4 — Neither alive 


] 
0 


is a matrix such that if z;is a column vector representing the relative prob- 
abilities at step 7, AX, will represent probabilities at step i+1. Then if 


we take X, = , X,, = A”X,. If we modify A as follows : 


0 
0, 
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‘ 
ps ] 0 0 
A= 
qr 0 0 
pr O ] 


A,=1, A,=1-6, A, =l-e, 


A, = 98 


then we can find eigenvectors which are solutions of: 


pr 


then we can find constants C; such that 


-) 0,0; 


r 8 
l-qs’ l-qs l-qgs 


(1-gs)(qs)” 


then 


1 ps-ps(qs)” 
1+gs\ gr—gr(qs)” 
1 - qs -gr-ps-pr(gs)” 


0 


8 8 

l-qs\ qr l-qs\ qr 
l-qs-gqr-ps pr 


The results can be expressed a little easier if we let 
K=l-qs, K, =1-(qs)" 
k( qs)” 


then 
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qs 0 0 
qr 1-5 0O 
pr € o 1 
- we can find eigenvalues which are solutions of |A-A/| = 0. : 
q 
Q. - 0 0. = 0 1 Ps 
1 0 v9 ] 0 v4 gr 
é = . 
Xn 
or 
lf psk,, x - lf ps 
K\ K\ 9 
rk, pr 
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Also solved by George M. Bergman, Stuyvesant High School, New 
York; Jeff Cheeger, Erasmus Hall High School, Brooklyn, New York; and 
Melvin Hochster, Stuyvesant High School, New York. Two incorrect solu- 
tions were received. 


A System of Equations 


380. [May 1959] Proposed by Huseyin Demir, Kandilli, Eregli, Kdz., Tur- 
key. 
Solve the system of equations 


(1) a(z—a)+ulr+u) =0 
(2) y(a— b)+uly+u) = 0 
(3) ay—c)+u2+u) =0 


1 1 


where abc # 0 and a~'+07'+c7'=u7?~ 

Solution by Chih-yi Wang, University of Minnesota. By performing the 
operations multiply (1) by y, multiply (2) by z, multiply (3) by 2 and ap- 
plying (2), (3), (1) respectively we get 


1 


ay2-aby+auy+au?+ buy-u® = 0 
= 0 
= 0 

By performing the operations (4)-(5), (5)-(6), (6)-(4) we get 

(7) (au + bu-ab)y+(be- bu-cu)z = (b-a)u? 

(8) (ca-cu-au)x+(cut+ bu-be)z = (c-b)u? 

(9) = (a-c)u? 


Since the augmented matrix of (7), (8), (9) is of rank 2, we can cal- 
culate two variables in terms of the third, so we get 


a) 
(10) 


b(c-5) 


ca 


(11) 


By substituting (10) into (3) we get, after simplification, 


0 
a 


whence by aid of (10) and (11), we obtain 
w2=(b/a)u, y=(e/b)u, 2=(a/c)u. 
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Note that we have used the relations abc #4 0, a~'+ b7'+e7! = u' when- 


ever necessary. 

Also solved by D. A. Breault, Sylvania Electric Products, Inc., Wal- 
tham, Massachusetts; Victor Ch’in, Kent State University, Kent, Ohio; 
Melvin Hochster, Stuyvesant High _~\ool, New York; and the proposer. 


Spheres About a Tetrahedron 


381. [May 1959] Proposed by George M. Bergman and Melvin Hochster, 
Stuyvesant High School, New York. 

Show that if the circumsphere and the insphere of a tetrahedron are 
concentric, four spheres can be drawn each tangent to each edge of the 
tetrahedron, extended if necessary. 

Solution by J. W. Clawson, Collegeville, Pennsylvania. Let A,, A,, 
A,, A, be the vertices of the tetrahedron; a,, the length of edge A,A,; 
B ,. the midpoint of that edge; O, the circumcenter of the face opposite 
A,; and so on; Othe circumcenter of the tetrahedron. 

Obviously OO , is perpendicular to the face opposite A, and 0,B,, to 
the edge A,A,; and so on, 

Now if O is also the in-center of the tetrahedron, 00, = 00,. Hence 


O = 9,B,,; triangle O ,.B,,A, = triangle 0,B,,A, and 0,A, = 0,A,. 
Thus the circumradii of the faces are all equal. It follows at once that face 


angles A,A,A,, A,A,A, are equal and in general that face angles oppo- 


site a common edge in two of the faces are equal. 

Consider triangle A ,A,A,. Let <A, = «, <A, = B, then <A, = 7-«-B. 
Then in triangle A ,A,A,, <A, = 7-«-f, let <A, = y, then <A, =«+B-y. 
In triangle A ,A,A,, <A, = B, <A, = y, then <A, = 7—B-y. And in triangle 
A,A,A,, <A, = 7-B-y, <A, = «+B-y, <A, = « Hence 7+ 2«-2y = 2, or 
« = y. It follows that all the faces are similar; and, since their circumradii 
are equal, that they are congruent. Thus opposite edges a,,, @,,; @,3, 


Q4° 14) Uy, are respectively equal. Note, however, the restriction that, 
since the sum of the face angles at any vertex is 7, and since the sum of 
two of these face angles must be greater than the third, all of the face an- 
gles are acute. Thus a tetrahedron whose circumcenter and incenter coin- 
cide has its opposite edges equal, but its congruent faces must be acute- 
angled triangles. 

Now, if a sphere touches the six edges of the tetrahedron, its sections 
by the four faces are circles inscribed or escribed to the triangular faces, 
Let the points of contact of the incircle to triangle A,A,A, be Pe jsp 


1,4, and those of the excircle opposite A, to triangle A,A,A, be ,J,,, 


1°38 
J J and so on, Let s be the semi-perimeter of each of the con- 


4% 28) 
gruent faces. 
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In 


Then Ag = 8 = = and A,,J,.= 8, Aa tog * 


this way we may show that the inscribed circle to any face touches any 
one of its sides at the same point at which the escribed circle to the other 
face containing that side touches that side; and that escribed circles to 
two faces which have an edge in common touch that edge produced in the 
same point. This proves the proposal. 

Also solved by the proposers. N. A. Court, University of Oklahoma 
submitted two bibliographical notes : 

a. If in a tetrahedron the circumcenter and the incenter coincide, each 
edge is equal to the opposite edge, and the tetrahedron is said to be isos- 
celes (Nathan Altshiller-Court, Modern Pure Solid Geometry, p. 97, art. 
304). 

b. With an isosceles tetrahedron four spheres may be associated each 
touching three edges internally and the remaining three edges externally 
(John Leech, Some Properties of the Isosceles Tetrahedron, Mathematical 
Gazette, London, vol. 34, no. 310 (December, 1950), p. 270, Theorem 8). 


The Hypocycloid 
382. [May 1959] Proposed by C. N. Mills, Sioux Falls College, South Da- 


ota. 
For the hypocycloid 7°“? +y?/% = a?/, («, 8) is the center of curva- 
1/8) 3, 


ture. Show that «+B = 
Solution by D. A. Breault, Sylvania Electric Products Inc., Waltham, 


Massachusetts. Since we have, 


(1) y = (a2/3 22/8) 


we find from a book of tables and formulas that the center of curvature, 
(«, 8), of y = f(z) is given by, 


(2) «= 2-y(l+y’?)/y”, 
(3) B 


whence 


(4) 


Calculation of y’ and y’’ from (1), and substitution into (4) gives the de- 
sired result: 


(5) «+B = 


(Intermediate calculations have been omitted to avoid making the proof 
very much longer than it is now.) 

Also solved by Norman Anning, Alhambra, California; Jeff Cheeger, 
Erasmus Hall High School, Brooklyn, New York; Victor Ch’in, Kent State 
University, Kent, Ohio; Melvin Hochster, Stuyvesant High School, New York; 
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F.D. Parker, University of Alaska; William M. Sanders, Mississippi South- 
ern College; C. M. Sidlo, Framingham, Massachusetts; Chih-yi Wang, Univ- 
ersity of Minnesota; Ron Wilder, Grants Pass, Oregon; Dale Woods, State 
Teachers College, Kirksville, Missouri, and the proposers. 


QUICKIES 


From time to time this department will publish problems which may be solved by la- 
borious methods, but which with the proper insight may be disposed of with dispatch. Read- 
ers are urged to submit their favorite problems of this type, together with the elegant sol- 
ution and the source, if known, 


Q 263. Solve arctan(p/z) + arc tan(g/x)+arctan(r/z) = 7 for where p, 
and r are fixed. [Submitted by Melvin Hochster and Jeff Cheeger| 


Q 264. Evaluate y > min(m,n) for M > N. [Submitted by M. S. Klamkin] 
m=1 N= 


Q 265. Show that S 1/p diverges. [Submitted by David L. Silverman] 
p prime 


Answers 
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TRICKIES 


A trickie is a problem whose solution depends upon the perception of the key word, 
phrase or idea rather than upon a mathematical routine. Send us your favorite trickies. 


T 34. Twenty sheep are to be slaughtered in five days such that (1) only 
an odd number of sheep can be slaughtered in each day, and (2) nokilling 
in any one day is ruled out since zero is naturally regarded as an even 
number. How can the schedule to slaughter be arranged? [Submitted by 
Chih-yi Wang] 

T 35. Double the area of a square hole in a wall without increasing the 
height or width of the hole. [Submitted by Lindley J. Burton] 


Solutions 
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Provides an introduction to the fundamental concepts 


and basic techniques of modern mathematical thinking 


Basic Concepts 
of Elementary Mathematics 


By William L. Schaaf, Brooklyn College. 


This work affords the reader a genuine insight into the precise nature 
of mathematics as an intellectual discipline. Rather than show how to per 
form mathematical operations, the author has chosen to show why these 
operations work; why multiplication, for instance, produces the answers it 
does, and why these answers are “correct”. This is done by presenting the 
unifying concepts and logical bases of the subject, and by demonstrating 
- how these are welded together into a rigorous deductive system. 

This whole approach follows from the notion that mathematics is not 
just a mass of Computational techniques but rather a way of thinking which 
can only be properly understood (and effectively taught) through an ac- 
quaintance with the basic assumptions of mathematical thought and the 
rules which govern it. With its lucid, informal, yet forthright style and its 
clarity of exposition, Professor Schaaf’s book makes the subject accessible 
to everyone interested in it. The book contains over 650 problems and ques- 
tions for the student to solve, as well as many references for further read- 
ing and study. 


1960 Approx. 380 pages Illus. Prob. $4.95 


Send for examination copies. 


JOHN WILEY & SONS, Inc. 
440 Park Avenue South New York 16, N.Y. 
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